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ABSTRACT 

The fundamental rovibrational band of CO near 4.7 /im is a sensitive tracer of the presence and loca¬ 
tion of molecular gas in the planet-forming region of protoplanetary disks at 0.01-10 AU. We present 
a new analysis of a high-resolution spectral survey (R^ 96,000, or ^ 3.2 km s -1 ) of CO rovibrational 
lines from protoplanetary disks spanning a wide range of stellar masses and of evolutionary properties. 
We find that the CO emission originates in two distinct velocity components. Line widths of both 
components correlate strongly with disk inclination, as expected for gas in Keplerian rotation. By 
measuring the line flux ratios between vibrational transitions F v= 2 -i/F v= i_q, we find that the two ve¬ 
locity components are clearly distinct in excitation. The broad component (FWHM = 50 —200kms _1 ) 
probes the disk region near the magnetospheric accretion radius at « 0.05 AU, where the gas is hot 
(800 — 1500K). The narrow component (FWHM = 10 — 50kms -1 ) probes the disk at larger radii 
of 0.1-10 AU, where the gas is typically colder (200-700 K). CO excitation temperatures and orbital 
radii define an empirical temperature-radius relation as a power law with index —0.3 ±0.1 between 
0.05-3 AU. The broad CO component, co-spatial with the observed orbital distribution of hot Jupiters, 
is rarely detected in transitional and Herbig Ae disks, providing evidence for an early dissipation of 
the innermost disk. An inversion in the temperature profile beyond 3 AU is interpreted as a tracer of 
a regime dominated by UV pumping in largely devoid inner disks, and may be a signature of the last 
stage before the disk enters the gas-poor debris phase. 

Subject headings: circumstellar matter — extrasolar planets — planets and satellites: formation — 
protoplanet ary disks — stars: pre-main sequence 


1. INTRODUCTION 

The inner 10 AU of protoplanetary disks is thou ght to 
be_the region in which most exoplanets form (Armitage 
2011). Initially, the inner disk likely has a smooth ra¬ 
dial distribution of gas and dust that is truncated by the 
stellar magnetic field at the co-rotation radius, within 
0.1 AU. From this point, disk material accretes onto the 
central sta r through magnetospheric accretion (e.g. Bou- 
vier et al.||2007|). However, over time, inner disks evolve, 
and g apsand inner holes fr equently develop within a few 
Myr (Muzerolle et al. 2010). S ome of these gaps may 
be carved by photo-evaporation (Alexander et al. 2014b, 
and others by dynamical clearing by planets coupled with 
pressure -driven dust migration (Zhu et al. 2011 Pinilla et 
al.||2015|). The evolutionary window available for planet 
Tormation is therefore limited by, and link ed to, the time 
scale of inner disk dispersal of a few Myr (Luhman et al. 
2010). When observing the inner regions of protoplane- 


tary disks and the process of planet formation, a central 
challenge is that the angular size of these regions is small 
(< 0.1"), even for the closest star-forming clouds. Fur¬ 
ther, the gas temperatures in the planet-forming region 
typically spans 100 — 1000 K, pushing the main obser¬ 
vational line tracers into the near- to far-infrared wave¬ 
length range. This leads to the common use of infrared 
molecular spectroscopy to investigate the properties of 
inner disk gas. 

The CO molecule is a ubiquitous and sensitive tracer of 
gas in protoplanet ary disks. While pure rotational CO 
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lines observed by millimeter interferometers trace cold 
gas extending to the outermost disk regions, CO lines 
from rovibrational transitions at infrared wavelengths are 
typically excited closer to the central star, where the disk 
is hotter and exposed to strong infrared (IR) and ul¬ 
traviolet (UV) radiation fields. Surveys of rovibrational 
CO emission from the fundamental (Av = 1) band near 
4.7 jam have demonstrated that the lines probe the prop¬ 
erties of inner disks between the mag netosp heric accre¬ 
tion radius and the inner 10 AU (e.g.jNajita et al.|[2 003 
Blake &; Boogert||2004| |Brittain et al.||2007| |Salyk et al. 


201ip. CO line profiles, when spatially- and velocity- 
resolved, have shown a wealth of CO isotopologues and 
complex line profiles, as indicative of a range in excita¬ 
tion cond itions and of physically distinct emission com- 
ponents (iPontoppidan et al. 2011b; Herczeg et al.||2011 

~ - al.|2oll| IJjrowiT ^ltrilM#" - ' - - N 1 


Bast et al.||2Ull| |Brown et al.||2013 |van der Bias et al 

2015p . Emission Rom —r7 UOu= 1-0 and v=2-l lines is 


strong enough to be detected in the majority of gas-rich 
disks, demonstrating that CO is often vibrationally ex¬ 
cited. In the innermost disk, CO can be efficiently ex¬ 
cited by IR photons from the stellar radiation and t he lo¬ 
cal /a m-size hot dust through IR pumping (e.g. Thi et al. 
2013). Detection of higher vibrational lines in some disks 
implies vibrational temperatures that are several times 
larger than the rotational temperatures, providing evi¬ 
dence for the radiative cascade from pumping of directly 
UV-irradiated CO gas (Krotko v et al.| 1980). Detection 
of rarer isotopologues, including 13 CO that is commonly 
seen in its v=l-0 transitions, has been usually attributed 
to large column densities of CO gas found in young disks. 

Recent analyses of velocity-resolved CO rovibrational 
emission in one protoplanetary disk found that the 
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Figure 1 . Portion of one of the 4.7/xm CRIRES spectra in our sample. Rovibrational P- and R-branch lines of 12 CO and 13 CO are 
marked in different colors. 


v— 1-0 lines are a superposition of a broad (FWHM 
150 km s -1 ) and a narrow (FWHM ~ 30kins -1 ) com- 


2015D- Her- 

y broad and 


pon ents (Goto et al. 2011; Banzatti et al. 
czeg et al. [i f201l| ) found evidence for similar 
narrow GO components in ^ 10 more embedded disks, 
and found that, while the two components are blended 
CO v=l-0 lines, only the broad component con- 


m 


12 


tributes to the 12 CO v=2-l lines, while the 13 CO lines 
are dominated by the narrow component. In this paper, 
we systematically apply a spectral decomposition analy¬ 
sis to high resolution 3 kms -1 ) spectra of a sample 
of ^ 40 protoplanetary disks, to test the hypothesis that 
CO rovibrational emission in protoplanet ary disks may 
generally be described by two distinct components in ve¬ 
locity and excitation: specifically, a narrow component , 
or NC, and a distinctly broad component , or BC. The de¬ 
composition allows us to estimate both the location and 
excitation temperature of the emitting CO gas. While 
most disks in the sample do show two distinct CO compo¬ 
nents that get colder at larger disk radii, some show only 
narrow emission lines and signatures of having cleared- 
out inner disk regions. Together, this dataset reveals a 
sequence of inner disk dispersal across a wide range of 
stellar masses, from 0.3 to 3 M 0 . 

2. OBSERVATIONS 

The disks sample includes 37 disks with 4.7 /am spec¬ 
tra obtained with t he C Ryogenic Infrared Echelle Spec¬ 
trometer (CRIRES, Kaeufl et al. 2004) on the Very Large 
Telescope (VLT) of the European Southern Observatory 
(ESQ), as part of the ES Q Large Program 179.C-0151 
Pontoppidan et al. (2011a). Further, one spectrum was 
added from program 093. C-0432, providing an additional 
epoch of data for a disk around the variable star EXLup 
(Banzatti et al. 2015). The original survey included disks 
mostly arouncTybung solar-mass stars (M* <1.5 M 0 , T 
Tauri stars) as well as some around intermediate-mass 
stars (1.5 < M* < 3 M e , mostly Herbig Ae/Be stars). 
The disks were also selected to cover a range of evolu¬ 
tionary phases, from the early stages when the disks are 
still embedded in diffuse envelopes to the clearing stages 
of transitional disks. For this paper, we selected those 
disks with spectra that show CO lines in emission, in¬ 
dicative of the presence of warm CO gas in the inner part 
of the disk. The sample includes 30 T Tauri disks and 7 
Herbig disks; 2 disks have been observed to have diffuse 
envelopes, and 7 disks are known to be transitional (see 
Table gj. 

All spectra were reduced using procedures described in 


Pontoppidan et al. (2011b), and were already presented 
in several papers, providing general overviews as well as 


more detailed studies of various snbsamples JPontoppi- 
dan et al.||20081 |2011bl pier czeg et al. ][20TTJ [B ast et al.^ 

. | Brown et al.||2013|j . The observations were carried 

between 2007 and 2010, but consistent instrumen- 


20TTT B 


out Detween 

tal parameters were maintained to ensure homogeneity 
of the data. The slit width was ~ 0.2”, providing a spec¬ 
tral resolution of ^ 3.2kms -1 , which is generally con¬ 
firmed by measured widths of narrow atmospheric ab¬ 
sorption lines. Early-type telluric standard stars were 
observed at similar airmasses to the science targets (typ¬ 
ically within 0.1 difference) to allow an accurate telluric 
correction. The observations were generally obtained us¬ 
ing the CRIRES adaptive optics correction. 

The four CRIRES detectors provide a combined in¬ 
stantaneous spectral coverage of ~ 0.1 /am. The observ¬ 
ing strategy targeted the 12 CO rovibrational P-branch 
lines in several non-contiguous settings, with at least 
one centered between 4.65 and 4.77 /am (see Figure [I]). 
Additionally, a spectral setting covering wavelengths up 
to 5.0 /im was obtained for most targets. The CRIRES 
spectra do not provide continuous coverage due to gaps 
from detector chips and the telluric absorption from at¬ 
mospheric CO lines Doppler-shifted from the targeted 
lines depending on the month of observation. In some 
cases, the spectra were repeated by applying a small shift 
to the central wavelength or taken in different months of 
the year to compensate for those eff ects and p rovide a 


more c omplete spectral sampling (see Pontoppidan et al. 


2011b). This strategy was adopted to obtain complete 
velocity-resolved CO line profiles. 


3. ANALYSIS OF CO SPECTRA 

In this section we describe the methodology applied 
to the analysis of CRIRES spectra, while comparison of 
our results to previous work and discussion are provided 
in Section [4] We apply a spectral decomposition algo¬ 
rithm, which essentially consists of the following three 
basic steps. 

Step 1 : High signal-to-noise profiles of 12 CO v=l-0 and 
v=2-l emission are built by stacking multiple lines 
from different rotational levels. 


Step 2: The stacked line profiles are split into two 
(broad; BC and narrow; NC) components by scal¬ 
ing the 12 CO v=2-l profile to match the wings of 
the 12 CO t>=l-0 profile, and then subtracting it. 
This step reveals if the emission is made of a single 
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Figure 2. Left: Empirical decomposition of CO rovibrational emission from protoplanetary disks into two components, where one (BC) 
is distinctly broader than the other one (NC). Comparison of 12 CO v=l-0 and r>=2-l lines shows whether two (double-component disks) 
rather than one (single-component disks) component is found in a given disk. Decomposition plots for the entire disk sample are provided 
in Appendix [b] Right: Comparison between 12 CO and 13 CO line widths. The 1:1 relation is marked with a solid line. 


(SC) rather than of two components (BC and NC, 
Figure [5]). 

Step 3: The line widths and the vibrational 
F v= 2 —\/F v= i—o ratios are then measured for 
each detected component. 


This procedure is applied to each CO spectrum in our 
sam ple. The first two steps are described in Sections |3.1| 
and |3.2| The third step provides the data used to find 
two empirical r elati ons that probe the properties of inner 
disks (Sections 3.3 and 3.5). 


3.1. Stacking of CO line profiles 

We build high signal-to-noise CO line profiles for each 
science target by stacking several transitions observed in 
the 4.7 /rm spectra. Nominally, the maximum spectral 
coverage obtained in this survey includes ro-vibrational 
CO transitions from P(l)-P(32) for 12 CO v— 1-0, P(l)- 
P(26) and R(0)-R(6) for 12 CO v=2-l, and P(l)-P(23) 
and R(0)-R(15) for 13 CO v=l-0. In practice, the num¬ 
ber of transitions available for each disk depends on the 
spectral coverage obtained for individual targets and on 
the gaps caused by detector chips and telluric absorp¬ 
tion, which in turn depends on the time of year of a 
given observation. For stacking and averaging line pro¬ 
files, we use the lines that are most commonly available, 
which have similar integrated flux, and which are not 
significantly blended with other transitions. These crite¬ 
ria preserve both the line flux and shape in the stacked 
profiles, and enforce homogeneity in the analysis method 
for the entire disk sample. Note that, within the rela¬ 
tively narrow spectral range 4.7 — 4.9 /im, lines out of a 
given vibrational level have very similar fluxes and exci¬ 
tation energies (within 10-30%) and line profiles depend 


only very weakly 

on the rotational quantum number (e.g. 

Salyk et al. 20llj 

Pontoppidan et al.|2011b;|van der Plasj 

et al. 2015). For 

high rotational quantum numbers, this 

; true (J > 25). 

may no longer be 


Given these constraints, the 12 CO 17=1-0 lines that are 
most suitable for stacking are those between P(7) and 
P(12), while the 12 CO v=2-l lines are those between P(3) 
and P(6). As for the 13 CO lines, several can typically be 


used over the entire spectrum, where not blended with 
12 CO. Each line is normalized to its local continuum, re¬ 
sampled on a common velocity grid, and the stacked line 
profiles are built by weighted average of the pixel values 
from the individual transitions. While the stacked v=2-l 
line profiles typically have continuous spectral sampling, 
v=l-0 lines often have spectral gaps on one side due to 
the counterpart telluric absorption lines at the Doppler 
velocity of the target at the time of observation (see for 
instance the gap on the red side of the ^=1-0 line in 
Figure J2|. In some disks (CW Tau, DoAr 24E S, HT 
Lup, VSSG 1), an absorption component is present on 
top of 12 CO emission in the t>=l-0 lines. This is typ- 
ically due to absorption from unrelated foreground gas 


iHerczeg et al.|2011 Pontoppidan et al.|2011b Brown et 

il.||2013|). In these targets, we exclude pixels affected by 


the absorption component from the stacked profiles and 
further analysis. 


3.2. Empirical decomposition into two CO components 

The presence of separate broad and narrow compo¬ 
nents is apparent by visual inspection of many disk spec¬ 
tra (e.g. Figure [I]). In practice, the difference is formal¬ 
ized by comparison of the t7=l-0 and v=2-l stacked line 
profiles. We scale the ^=1-0 to match the broad wings of 
the v=2-l profile, and then subtract one from the other. 
The best scaling factor is found by y 2 minimization be¬ 
tween the scaled v=l-0 line profile and the v=2-l line 
profile, excluding the central region of narrow emission, 
if present. The NC is typically much narrower than BC, 
such that most of the line profile can be used for the y 2 
minimization. The choice of the wavelength exclusion re¬ 
gion is determined iteratively, starting by using the full 
line profile. If the residuals after subtraction of the two 
line profiles are consistent with zero over the entire spec¬ 
tral range of the v=2-l line, the NC is not detected and 
the procedure stops there. In the contrary case, the fit 
is iterated after exclusion of increasingly large portions 
at the center of the line profile, until the residuals are 
composed of a flat continuum and of the NC emission 
line. Figure |2] shows a demonstrative sketch of the out¬ 
come, while spectral decomposition plots for all disks are 
shown in Appendix |B) 
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Table 1 

Sample properties and disk inclinations. 


Name 

CO a 

Type b 

M* 

[Mq] 

Ref 

j. c 
t, m 

[°] 

Ref 

iout c 
[°] 

Ref 

Comment 

AATau 

S 

TT 

0.85 

ml 

_ 

_ 

71(1) 

il 


AS205 N 

D 

TT 

1.1 

m2 

20(5) 

i2 

25 

i3 


AS209 

D 

TT 

1.4 

m3 

- 

- 

38-40 

i3 


CrA-IRS2 

D 

TT,Em 

1.0 

m4 

17(11) 

this work 

- 

- 


CVCha 

D 

TT 

2.1 

m3 

30(12) 

this work 

- 

- 


CWTau 

D 

TT 

1.2 

m2 

28(10) 

this work 

- 

- 

Zi n agrees with Zj et =41 from i4 

DFTau 

D 

TT 

0.5 

m5 

65(15) 

this work 

- 

- 

ii n agrees with z s tar = 60 from i5 

DoAr24E S 

D 

TT 

0.7 

m2 

20(5) 

i2 

- 

- 


DoAr44 

S 

TT,Tr 

1.3 

m3 

- 

- 

40(20) 

i3 


DRTau 

D 

TT 

1.0 

m2 

9(5) 

i2 

37(3) 

i6 


EC82 

S 

TT 

1.0 

m4 

51(18) 

this work 

- 

- 


EXLup 08 

D 

TT 

0.6 

mb 

44(8) 

this work 

- 

- 

EX Lup observation from August 2008 

EXLup 14 

D 

TT 

0.6 

m6 

48(17) 

this work 

- 

- 

EX kup observation from April 2014 

FNTau 

S 

TT 

0.35 

m7 

- 

- 

20(10) 

i7 


FZTau 

S 

TT 

0.7 

m3 

38(15) 

this work 

- 

- 


GQkup 

D 

TT 

0.8 

m2 

65(10) 

i2 

- 



HD135344B 

S 

H,Tr 

1.6 

m2 

14(4) 

i2 

23(0.3) 

i8 


HD142527 

S 

H,Tr 

3.5 

m2 

20(2) 

i2 

23(8) 

i9 


HD144432S 

s 

H 

1.7 

m2 

27(3) 

i2 

- 

- 


HH100 

D 

TT,Em 

0.4 

m3 

26(4) 

this work 

- 

- 


HTLup 

D 

TT 

2.5 

m3 

28(8) 

this work 

- 

- 


IMLup 

D 

TT 

1.0 

m4 

- 

- 

49(4) 

ilO 


IRS48 

S 

H,Tr 

2.0 

m3 

42(6) 

ill 

48(8) 

il2 


LkHa330 

s 

H,Tr 

2.5 

m2 

12(2) 

i2 

42 

il3 


RNO90 

D 

TT 

1.5 

m2 

37(4) 

i2 

- 

- 


RULup 

D 

TT 

0.7 

m2 

35(5) 

i2 

- 

- 


SCrA S 

D 

TT 

0.6 

m3 

26(4) 

this work 

- 

- 


SCrA N 

D 

TT 

1.5 

m2 

10(5) 

i2 

- 

- 


SR21 

S 

H,Tr 

2.2 

m2 

15(1) 

i2 

15(0.4) 

i8 


TTau N 

D 

TT 

2.1 

m3 

20(13) 

this work 

- 

- 

Zi n agrees with Zdisk < 30 from H4 

TWCha 

S 

TT 

1.0 

m8 

67(23) 

this work 

- 

- 


TWHya 

S 

TT,Tr 

0.7 

m2 

4(1) 

i2 

12 

i6 


VSSG1 

D 

TT 

0.5 

m3 


- 

53 

i3 


VVSer 

S 

H 

3.0 

m2 

65(5) 

i2 

70(5) 

H5 


VWCha 

D 

TT 

0.6 

m3 

44(17) 

this work 


- 


VZCha 

D 

TT 

0.8 

m3 

25(10) 

this work 

- 

- 


WaOph6 

D 

TT 

0.9 

m9 

- 


39 

i3 


WXCha 

D 

TT 

0.6 

m8 

87(31) 

this work 

- 

- 



References. — Stellar masses: (ml ) |Bouvier et al.| (|1999j), (m2) |Pontoppidan et al.| (|201 lbj) , (m3 ) [Salyk et 
ak| ([2013]) an d refere nces therein, (m4) assumed eq ual to t he median stellar mass in this sample, (m5) |B< 
et al.N l99U|, (m6) Gras-Velazquez Ray ( 2005| ), ( m7) [ Mc Clure e t al. (2013), (m8) Feigelson et al.' 

t r\4- n 1 \ or\L i /-» c\ 1 n r\ 


eckwith 


(my) Andrews et 

Pontoppidan et al. |~f 

%4 ) Coffey et al. 


Disk inclinations: 


Cox et al." 


image, 


|Q8 | |201 lb|) from CO spectro-astr omet r. 

[2008p , (z5) Unruh et al.| ( |1998|, (i6 ) Isella et al. (2009) 


2013) ~Irom optical scatter ed light image, (z 
Andrews et al.| (|2009|) from resolved mm 
from resoiv 


mid-infrared (MIR) image 
resolved MIR image, (i 


et al.| (|zui4|) irom i Nm s cattered lignt image, \i±U) rmte et al. (|z(JU8|) trom resol ved mm and op 
light images, (ill) Brown et al. ( 2012|) from r esolve d 1N1R CO image, (H2) [Geers et ak] (|2007|) 


ilved mm image, (z7) 


Kudo et ak (2008) from KIR scattered light image, (fe] Perez et ak (2014)mom resolved mm image, (zP)|Avenhaus 
et al.| ( |2014| from NIR scattered light i mage, (HO) Pint e et al.j (J 2008) from resolved mm and op tical scattered 


age, [il3) Brown e t al.| (|2009|) 
15) |PontoppI dan et al.f( 2007f : 


_ , , from resolved 

from resolved mm image, (il4) Ratzka et ak (2009) from 
from resolved MIR image. 


Note. — 

a Indication whether two (D) or one (S) component of CO emission are detected. 

b T Tauri stars (TT) have M* < 1.5 Mq, and Herbig Ae/Be stars have M* > 1.5M 0 . Disks embedded in diffuse 
envelopes (Em) or defined as transitional disks by previous work (Tr) are marked accordingly. 
c Disk inclinations are reported for both the inner and the outer disk, depending on the technique used. Errors 
are provided in brackets, where available. 


Two basic parameters are measured from the decompo¬ 
sition procedure: the line width and the vibrational flux 
ratio v2/vl = F v= 2 -i/F v= i_q. Line widths are mea¬ 
sured at 50% of the line peak (the full width at half 
maximum; FWHM) from fits to the line profiles, allow¬ 
ing up to two Gaussian functions to account for double- 
peaked profiles. The vibrational flux ratio is measured 
from the scaling factor that matches the 17=1-0 to the 
v—2-1 lines as explained above, and is exactly equivalent 
to the integrated line intensity ratio between the v=2- 


1 and the v=l-0 lines in each spectrum (further details 
are provided in the Appendix). It is possible to mea¬ 
sure vibrational ratios in all single-component disks and 
for the BC in all double-component disks. For the NC, 
F v = 2 -i/F v= i-o can be measured in only three disks: AS 
205 N, DR Tau, and S CrA N. In these disks, the NC 
can also be identified in the v=2-l lines through fits of 
multiple Gaussians to the line profile, one for each of the 
broad and narrow components (plots are shown in the 
Appendix). The non-detection of the NC in the v=2- 
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Table 2 

Disks with two CO components. 


Table 3 

Disks with a single CO component. 


Name 

BC 

FWHM 

[kms -1 ] 

v2/vl 

Rco 

[AU] 

^vib.ex. 

[K] 

NC 

FWHM 

[kms" 1 ] 

Rco 

[AU] 

AS205 N 

52 

0.38 

0.17 

970 

13 

2.3 

AS209 

160 

0.91 

0.07 

1054 

33 

1.7 

CrA-IRS2 

53 

0.37 

0.11 

960 

17 

1.1 

CVCha 

102 

0.25 

0.18 

790 

47 

0.9 

CWTau 

54 

0.32 

0.32 

890 

45 

0.5 

DFTau 

89 

0.43 

0.20 

1060 

53 

0.6 

DoAr24E S 

81 

0.30 

0.04 

860 

29 

0.4 

DRTau 

37 

0.39 

0.06 

990 

13 

0.5 

EXLup 08 

155 

0.61 

0.06 

1420 

32 

1.3 

EXLup 14 

150 

0.51 

0.07 

1220 

40 

1.0 

GQLup 

107 

0.28 

0.20 

830 

62 

0.6 

HH100 

80 

1.02 

0.04 

1200 

12 

1.9 

HTLup 

190 

0.65 

0.05 

1510 

35 

1.6 

IMLup 

200 

1.16 

0.05 

1400 

40 

1.3 

RNO90 

92 

0.21 

0.23 

730 

52 

0.7 

RULup 

111 

0.50 

0.07 

1200 

18 

2.5 

SCrA S 

101 

0.40 

0.04 

1000 

15 

1.8 

SCrA N 

57 

0.30 

0.05 

860 

12 

1.1 

TTau N 

110 

0.57 

0.07 

1340 

28 

1.1 

VSSG1 

66 

0.36 

0.27 

940 

35 

1.0 

VWCha 

80 

0.32 

0.16 

890 

33 

0.9 

VZCha 

53 

0.44 

0.18 

1080 

24 

0.9 

WaOph6 

133 

0.33 

0.07 

900 

32 

1.2 

WXCha 

140 

0.45 

0.11 

1100 

67 

0.5 


Note. — Rco is estimated from Kepler’s law using the line 
velocity at HWHM, the stellar mass, and the disk inclination 
(Table [TJ. XUb.ex. is estimated from v2/vl — F v =2-i/Tv=i-o 
assuming a CO column density of 10 18 cm -2 , apart from for 
AS209, IMLup, and HH100 where 10 19 cm -2 is used (see the 
Appendix for details). v2/vl and T v ib.ex. can be measured for 
NC only in three disks, where it is detected also in the 12 CO 
v=2-l lines: AS205 N (0.06, 520 K), DRTau (0.13, 620 K), and 
SCrA N (0.06, 520 K). 


1 lines of all other double-component disks is generally 
consistent with an upper limit on the vibrational flux 
ratio of < 0.2. Tables [2] and [3] report line widths and 
vibrational ratios for the entire disk sample. 

By implementing this decomposition procedure, we 
find that all disks fall into two phenomenological classes: 
the “single-component” disks (14 of 37), where the v—1- 
0 line profile, once scaled down, is consistent with that of 
the v=2-l line, and the “double-component” disks (23 of 
37), where the v=2-l profile matches only the wings of 
the 17=1-0 profile, revealing an additional narrower com¬ 
ponent of the t;=l-0 profile. In the double-component 
disks, the 13 CO lines are typically significantly narrower 
than the BC, but match the profile of the NC (Figure 
[2|. Only AS 209 and IM Lup have 13 CO lines de¬ 
tected in the BC only, rather than in the NC, while 
a few disks tentatively show both the BC and NC in 
13 CO (CrA-IRS 1, AS 205 N, and RU Lup). In double¬ 
component disks, the BC profile shows properties distinct 
from those of the single-component disks: broader lines 
(~ 50-200 km s -1 ), a clear separation in terms of line 
flux ratios with a boundary at F v= 2 -i/F v =i-o ~ 0.2, 
and weak or undetected 13 CO lines. Conversely, the NC 
shares the same line widths (~ 10-70 km s -1 ), same vi¬ 
brational ratios F v= 2 ~i/ F v= i-q < 0.2, and same 13 CO 
detection rate as for CO in 12/14 single-component disks 
(SC) of our sample. 


Name 

FWHM 

v2/vl 

Rco 

^vib.ex. 


[km s _1 ] 


[AU] 

[K] 

AATau 

115 

0.14 

0.2 

630 

DoAr44 

72 

0.11 

0.4 

590 

EC82 

49 

0.03 

0.9 

420 

FNTau 

9 

0.01 

1.8 

330 

FZTau 

31 

0.20 

1.0 

720 

HD135344B 

14 

0.05 

1.7 

420 

HD142527 

25 

0.04 

2.3 

460 

HD144432S 

40 

0.09 

0.8 

560 

LkHa330 

9 

0.09 

4.7 

420 

TWCha 

65 

0.16 

0.7 

660 

TWHya 

6 

0.09 

0.3 

500 

VVSer 

55 

0.16 

2.9 

660 

IRS48 

14 

0.38 

16. 

970 

SR21 

9 

0.64 

6.5 

1480 

HD100546 * 

16 

0.27 

15. 

810 

HD97048 * 

17 

0.33 

14. 

900 

HD179218 * 

18 

0.37 

21. 

950 

HD141569 * 

16 

0.56 

17. 

1300 

HD190073 * 

16 

0.27 

6.0 

810 

HD98922 * 

21 

0.29 

9.3 

840 


Note. — Rco is estimated from Kepler’s 
law using the line velocity at HWHM, the stel¬ 
lar mass, and the disk inclination (Table pi. 
Tvib.ex. is estimated from v2/vl assuming a CO 
column density of 10 18 c m~ 2 . Disks marked 
with * are included from Ivan der Plas et al.l 
(2015). 


3.3. An empirical relation between CO line width and 
disk inclination 

When using line widths to infer radial properties of 
Keplerian disks, it is necessary to correct for inclina¬ 
tion effects. Consequently, we collect values of disk 
inclinations from the literature, where available, for 
the disks in our sample. We carefully searched for 
the most reliable estimates obtained using spatially- 
resolved disk images, either from millimeter interferom¬ 
etry, optical/near-infrared (NIR) scattered-light imag¬ 
ing, or spectro-astrometry from CRIRES. Note that for 
spectro-astrometry the derived disk inclination is degen¬ 
erate with an assumed stellar mass. For consistency, we 
therefore adopt both the stellar masses and the disk in¬ 
clinations reported by|Pontoppidan et al.|(|2011b|). Table 
[l] shows the adopted disk inclinations, with notes on the 
technique used to measure the inclination for each disk. 
Using ten double-component disks with known inclina¬ 
tions (seven from spectro-astrometry, and three from in¬ 
terferometry), we find strong correlation between the CO 
line width and disk inclination, with Pearson coefficients 
of 0.7 for the BC and of 0.9 for the NC. Strong corre¬ 
lation between NIR CO line widt hs and disk inclination 
has been reported previously by|Blake & Boogert (2004) 
from a sample of five Herbig disks. Correlation between 
these parameters is expected when the gas is in Keplerian 
rotation around the star, a condition that is suggested by 
the high detection rate of double-peaked line profiles in 
our disk sample (> 70%) and that has been proposed 
since the first surveys of NIR CO emission in disks (e.g. 


Najita et al. 2003[ [Brittain et al. 


t 2003JL 

2004). Note that the two-component decomposition sig¬ 

nificantly decreases the scatter in this relation relative to 


Blake fc Boogert) 
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Known inclinations 



Derived inclinations 



Single-component disks 



0 20 40 60 80 

Disk incl [deg] 


Figure 3. Empirical relations between CO line widths and disk 
inclinations. Top: Linear fits to the distributions of the BC and 
the NC in disks with known inclinations from in the literature, 
based on spatially-resolved images. Middle : Best-fit disk incli¬ 
nations derived for disks with no direct inclination measurements 
from imaging, using the two linear relations from the plot above. 
Bottom : The single-component disks. Grey dashed curves indi¬ 
cate models for gas in Keplerian rotation at three different disk 
radii around a one solar-mass star, as indicated near each curve. 


what is seen when the line profiles are not fully velocity- 
resolved as with CRIRES, or if the line profiles were all 
modeled with a single Gaussian profile to derive average 
line widths (i.e. without separating BC from NC). 

Gas in Keplerian rotation at a disk radius R produces 
emission at a maximal velocity Av from the rest fre¬ 
quency of a given CO line following Kepler’s law as 


Av = sin(z) yjGM */R, where i is the disk inclination, G 
is the gravitational constant and M * is the stellar mass. 
Using this relation, it is possible to derive an emitting ra¬ 
dius for the gas in a rotating disk from the observed line 
widths. In general, emission from higher velocities orig¬ 
inate from smaller orbital radii. These velocities, usu¬ 
ally taken at the 50% or 10% of the line peak value, 
therefore provide estimates of the disk radii where CO 
is emitting. We refer to the innermost emitting radius 
of CO, Ri n , when derived from 1.7xHWHM as proposed 
by Salyk et al. (2011). In comparison, the velocity at 
l.OxHWHM provides an estimate of the disk radius Rco 
that emits the peak line flux (in agreement with t he em it¬ 
ting radii derived from spectro-astrometry by Pontopp- 
idan et al. 2008 2011b| ). Figure [3] shows the measured 


CO line widths divided by y/M*. Removing the stellar 
mass dependence on the inclination relation leaves only 
the dependence on the Keplerian radius R, which can be 
compared to a simple model of gas emission in Keplerian 
rotation around a 1 M 0 star (shown in Figure [3] at three 
different disk radii). 

In these terms, the BC has R in between 0.01 and 0.1 
AU. This is consistent with the disk region of the coro¬ 
tation and mag netosp heric accretion radii for solar-mass 
stars (e.g. Bouvier et al. 2007), as found by previo us stud¬ 


ies of rov ibrational CO emission from disks (e.g. Najita 


et al.|2003 Salyk et al. 2011). The NC, instead, emits 
from larger disk radii, with R in between 0.1 and 1 AU. 
The single-component disks show CO emission at radii 
similar to, or larger than, those of the NC. 

In Figure [3] we show linear fits to the BC and NC 
in disks with known inclinations, obtained using the 
Bayesian method by Kelly (2007). The regression takes 
into account the uncertainties on both variables, as well 
as the intrinsic scatter in the regression relation due to 
physical quantities not explicitly accounted for. Typical 
errors on the line widths are < 7kms -1 for the BC and 
< 4kms -1 for the NC, the stellar masses are assumed 
to be known to 20% and the disk inclinations are typi¬ 
cally known to accuracies better than 10%. The best-fit 
parameters and their errors are taken as the median and 
standard deviation of the posterior distributions of the 
regression parameters. The median intercept and slope 
are used to plot the lines in Figure [3J while the posterior 
distributions are plotted in the Appendix. 

The single linear fits shown in Figure [3] assume that 
differences in R in due to the dependence on stellar and 
disk luminosity for either the BC and the NC is weak 
relative to the inclination dependence. This is supported 
by the Bayesian regression of line widths versus disk in¬ 
clinations, which confirms strong correlation (> 0.8) for 
both the BC and NC, although the uncertainties in the 
regression parameters for BC are 3-4 times larger than 
for NC. This is also consistent with previous results for 
disks around solar-mass stars w here R\ n was fo und to 
vary only within a factor of a few (Salyk et al. 20lll), and 
with the theoretical expectation that the inner radius in¬ 
creases only with yfLl. The residual intrinsic scatter on 
the linear relation for the BC is four times larger than 
for the NC, which may point to a remaining dependence 
on the stellar and (variable) accretion luminosity. 
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CO vibrational ratio 



CO temperature 


3.2 ~ 


3.0 


> 2.8 - 


o 


2.6 
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log R co [AU] 


Figure 4. Left: Empirical relation found between vibrational flux ratios F v= 2 -i/F v= ±-o and the CO orbital radii Rco • The data 
include all CO components, broad or narrow, detected in the entire disk sample. Right: As the left panel, but with the vibrational flux 
ratio converted into an excitation temperature as explained in the text. This relation is essentially an empirical radial profile for the CO 
gas temperature in inner disks. SR 21 and IRS 48 are excluded from the fits. 


3.4. New estimates of disks inclinations 

Given the empirical line width-inclination relation, it 
is possible to estimate inclinations for disks with high 
resolution CO spectroscopy even when resolved imag¬ 
ing is not yet available. In doing so, it is assumed that 
the two relations found for the BC and NC represent 
a general property of rovibrational CO emission in the 
double-component disks. The advantage of having two 
relations, one for each velocity component, is that for a 
given double-component disk there are two constraints 
to derive one parameter. We take the measured values 
of the BC and NC line widths in each disk and use them 
to derive the disk inclination by using the regression co¬ 
efficients a, b from the empirical linear relations: 

i = (FWHM /y/Ml - a)/b, (1) 

where a = 31 ± 30 and b m 1.8 ± 0.9 for the BC, and 
a = 0.6 ±8.8 and b = 1.0 ±0.3 for the NC. We derive two 
inclination values, one for the BC and one for the NC, 
and adopt the mean of the 1 -a range in common between 
the two. Among the single-component disks, only three 
have unknown inclination: EC 82, FZ Tau, TW Cha. 
For these disks, we adopt the NC relation, since single¬ 
component disks with known inclinations span disk radii 
similar to those of the NC of double-component disks 
(Figure [3|. Following this procedure, we estimate an 
inclination for a total of fifteen new disks. These values 
are included in Table [l] 

3.5. A diagram of CO vibrational ratios and orbital 

radii 

Using the inclinations described in the previous sec¬ 
tion, we derive the CO emitting radii for all disks in 
the sample. In Figure [4] we show a diagram made with 
the vibrational line flux ratios plotted against the CO 
emitting radii. The diagram at the left of the Figure 
shows a strong logarithmic anti-correlation between Rco 
and F v= 2 -±/F v= i_q, with a Pearson coefficient of -0.9. 
In addition to spanning spatially separated Rco ranges, 
as already shown in Figure [3j the BC and NC clearly 
span distinct values in F v= 2 -i/F v= i_ 0 . Here too, the 
single-component disks have similar values to those of 


the NC. Two obvious outliers are SR 21 and IRS 48, 
both transitional disks with large inner gaps, which have 
vibrational excitation temperatures as high as the BC at 
much smaller disk radii (they will be discussed in Sec¬ 
tion [4|. Excluding them, Figure [4] demonstrates that the 
vibrational flux ratio decreases with disk radius as: 


F v =2-1 __ 7?—0.7±0.1 


F v = i-o 


oc R 


( 2 ) 


The vibrational line flux ratios measure how much the 
vibrational level v=2 is populated as compared to u=l, 
which in turn is a measure of the v ibrational excitation 
temperature T vi b.ex. (e.g. Thi et al.|[20l3 ). The relation 


between i^= 2 -i/^=i-o and Rco can thus be converted 
into a radial profile for the CO vibrational excitation 
temperature in inner disks. The conversion of the vi¬ 
brational flux ratio to a temperature depends on the 
CO column density, so the model-independent relation 
is that using the vibrational flux ratio (Eqn. 2). To es¬ 
timate the corresponding temperature relation, we use 
a model of a slab of gas in thermal equilibrium, where 
the level populations are defined by a single tempera¬ 
ture following the Boltzmann distribution. Other analy¬ 
ses of CO rovibrational emission from disks consistently 
find a relatively narrow range of column densities. Con¬ 
sequently, we adopt a single 12 CO column density of 
10 18 cm -2 , as an order-of-magnitude value representative 


±7 


± 9 , 


— 2 


Blake & Boogert 20041 IHerczeg et al. 2011 

Salyk et al. 

2011; Banzatti et al. 2015 van der Plas et a 

. 2(J15j). The 

details of the conversion are described in tJ 

ie Appendix 


|A| The resulting temperature radial profile is: 


T vib . ex . (X ir 0 - 3 * 0 - 1 (3) 

If disks without previously known inclinations are ex¬ 
cluded, the exponent becomes —0.2 ±0.1, within the er¬ 
ror, indicating that the result is not biased by our empir¬ 
ical derivation of disk inclinations. The excitation tem¬ 
peratures estimated for each disk are reported in Tables 
[2] and |3j The nominal uncertainty on the temperatures, 
once the column density is fixed, is less than 50 K as 
due to the uncertainty on the vibrational ratios (which 































Table 4 

Summary of rovibrational CO emission properties. 


CO comp. 

FWHM 

[kms -1 ] 

v2/vl 

Rco 

[AU] 

^vib.ex. 

[K] 

13 co 

Broad (BC) 

50-200 

0.2-0.6 

0.04-0.3 

800-1500 

3/23 

Narrow (NC) 

10-70 

< 0.2 

0.2-3 

< 700 

18/23 

Single I (SC) 

6-70 

0.01-0.2 

0.2-5 

300-700 

9/12 

Single II (SC) 

9-20 

0.3-0.6 

6-20 

800-1500 

8/8 


Note. — Single-component disks are separated into two 
groups based on their vibrational ratios, with F v = 2 -i / F v =i-o < 
0.2 found for the IR regime and 0.3 < F v = 2 -i/F v = p_o < 0.6 
found for the UV regime shown in Figure [5] including the 6 
Herbig disks from van der Plas et al. ( 2015| ). The last column 
reports 13 CO detection fractions. 


is typically < 10%). 

Herczeg et al. ( 2011 1 ) proposed that the non-detection 


of 13 CO lines in the BC, combined with the curvature 
seen in rotation diagrams, constrained the BC column 
density to within 10 17,5 -10 18 cm -2 . On the other hand, 
rotation diagrams of 13 CO lines were found to require 
column densities of 10 18,5 -10 19 cm -2 for the NC. If we 
relax our assumption and allow the BC to have columns 
as low as 10 17 cm -2 , and NC as large as 10 19 cm -2 , then 
the temperatures will decrease/increase for the NC/BC, 
respectively, and the overall temperature profile may ex¬ 
hibit a steeper index up to q ~ 0.5 (see Appendix |A| . 

4. DISCUSSION 

4.1. Origin of the broad and narrow CO components in 

disks 

The vibrational excitation temperatures of the BC 

match those previously found from rotational analyses of 
12 


CO lines (700-1700 K, e.g. 


Salyk et al.|2011 ), while the 


CO in the NC and SC disks typically have lower temper¬ 
atures, corresponding to those determined for 13 CO line 


fluxes (200-600 K, e.g. Brown et al. 2013). This is con¬ 


sistent with our finding that the 13 CO line profiles match 
those of the NC, as indicative of probing a similar disk 
region. In turn, this suggests that the NC is formed in an 
optically thick layer, leading to high F( 13 CO)/F( 12 CO) 
ratios, while the BC must be more optically thin. A sim¬ 
ilar c onclusion was r eached for disks in embedded young 
stars ( Herczeg et al.| 2011). 

The radii derived from the BC span the range found 
previously for disks around solar-mass stars (0.01- 
0.1 AU) while the NC and SC span the range found 
for Herbig and transitional disks (0.1-10 AU). This is 
consistent with the finding that double-component disks 
(where BC is detected) are primarily disks around low- 
mass stars, while the single-component disks are primar¬ 
ily Herbig and transitional disks. The NC and part of SC 
share similar properties in terms of line widths, orbital 
radii, and vibrational ratios (Section [3] and Table®, and 
they may trace the same gas component, regardless of 
the presence of a BC. 

While a disk origin for the BC is supported by the com¬ 
mon double-peak profiles, as well as b y a Kepl eri an signa l 
seen in spectro-astrometry (Pontoppidan et al.||2011b), 
the origin of the NC is less clear. A disk wind has been 
proposed for single-peak 12 CO i?=l-0 line profiles (those 
dominated by the NC, e.g. AS205 N) to reproduce the 


sub-Keplerian motio ns ofM he line p eak, o bserv ed wit h 
spectro -astrometry (|Bast et al.||2011| |Pontoppidan et al. 
2011b). In one case, a corresponding lar ge-scale disk 
wind nas possibly been seen with ALMA (Salyk et al. 
2014). An observational characteristic of a disk wind 


launched at small radii, and preserving angular momen¬ 
tum as it expands, is that it the gas contributing to the 
line becomes increasingly sub-Keplerian closer to the line 
center. Specifically, the line wings are exactly Keplerian, 
while the inner few kms -1 are strongly sub-Keplerian, 
tracing radii significantly smaller than what their veloc¬ 
ity would suggest. This could explain why observables 
constructed from the line wings of both the broad and 
narrow components, as we do here, yield results consis¬ 
tent with a nearly Keplerian disk. 

Broad and narrow components that may be similar to 
those of the rovibrati onal CO lines have be en found in 
the 6300A[OI] line by Rigliaco et al. (2013). These au¬ 
thors suggest a disk origin for the broad component, and 
an unbound cool disk wind as the origin of the narrow 
component. It is not clear whether the narrow compo¬ 
nents of [OI] and CO are generally tracing the same disk 
wind, but they share close resemblance in at least a few 
disks. For the disks in common between this work and 


Rigliaco et al. (2013), the FWHM of NC in CO and [OI] 
agrees to within < 3kms -1 in DR Tau, VZ Cha, RU 
Lup, TW Hya, and CW Tau; in DR Tau, the same is 
true also for the BC of CO and [OI]. 

4.2. The radial temperature profile as a tracer of 
excitation mechanisms 


In Section 3.5 we found a strong logarithmic anti¬ 
correlation between the CO vibrational excitation tem¬ 
perature and emitting radius. This relation corresponds 
to an empirical temperature profile of a typical inner disk 
surface. However, it was noted that two disks from the 
sample with large CO emitting radii, the transitional 
disks SR 21 and IRS 48, clearly depart from the rela¬ 
tion, showing significantly higher vibrational ratios at 
their respective CO emitting radii than the relation pre¬ 
dicts. To explore this further, we add 6 disks from the 
CRIRES survey of 12 disks around intermediate-mass 
2 — 3 Mq) stars from 


der Plas et al. 


(2015) 


(M* 

This survey presents 4.7 /am spectroscopy with a simi 
lar experimental design and wavelength coverage as our 
survey. Specifically, we include the full subset of disks 
for which lines are detected from both 12 CO v=l-0 and 
v=2-l transitions in a wavelength range similar to that 
of our survey. Following the same procedure as for other 
disks, we derive Rco and T v ib e x for the disks around 
HD100546, HD179218, HD190073, HD97048, HD141569, 
HD98922 (included in Table [^J). Their locations in the 
temperature-radius diagram are included in Figure [5| 
which shows that disks around intermediate-mass stars 
have vibrational temperatures that are typically higher 
at a given radius than the relation defined by the lower- 
mass stars at smaller disk radii. It is clear that SR 21 
and IRS 48 belong to this vibrationally hot class of disks. 
The combination of CRIRES surveys of disks spanning a 
wide range of stellar masses therefore demonstrates the 
existence of two trends: a vibrational temperature de¬ 
crease out to radii of ~ 3 AU followed by a temperature 
increase at radii > 3 AU. We note that the remaining six 
disks from van der Plas et al. (2015) lack v=2-l mea- 
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10 


R co [ AlJ ] 


Figure 5. CO temperature-radius diagram constructed by com¬ 
bination of CO rovibrational spectra from this work and spectra 


from the survey by | van der Plas et al. ( |2015jt The datapoints from 


this work are the same as shown in F 


surements due either to poor telluric corrections and low 
signal-to-noise, or the line widths and disk inclinations 
put them in the IR regime of our diagram (< 3AU), 
where low n=2-d^lrie_jAixes_are ; expected^ 


Ana lyses by Brittaine^ah] (2003), Blake fc Boogert] 


(2004), and van der Bias et al. (2015) included two pos¬ 


sible excitation regimes for CO rovibrational emission: 
in the inner disk, IR pumping due to the local dust 
continuum radiation, and at larger disk radii, where 
the dust temperature is too low to generate sufficient 
amounts of 4.7 /im photons, pumping by UV photons 
from the star and/or accretion shock. The specificity of 
UV pumping is to populate very high vibratio nal levels 


(v > 6) by decay of an excited electronic state (Krotkov 
et al. 1980), while IR pumping drives the CO level pop¬ 


ulation to equilibrate to the local color temperature of 
the dust. At dust temperatures of 500-1000 K, this will 
only populate the v=l and v=2 levels. Indeed, Figure [5] 
can be interpreted as evidence for a transition between 
the IR and UV pumping regimes near 3AU. The “IR 
regime”, at small disk radii, and the “UV regime”, at 
larger disk radii, are correspondingly marked in the 
Figure. Further evidence is provided by comparison 
to rotational temperatures T rot .ex . derived in the same 


disks from previo us studies (Herczeg et al. 
et al.||2QQ9 , 2011| ). Figure [5] shows that 'lAb, 


2011 


Salyk 


in the IR regime, indicating that the CO rovibrational 
population is ther malized to the local dus t temperature 
in the inner disk (Blake & Boogert 20041). Conversely, 


in the UV regime, 1 i-ot.ex. < -^vib.ex. (jThi et al. 
This is consistent with the finding of others, 
stance 


2013) 


Brittain et al. 


or m- 
= 200 K 


(2003) found T rot . ex , 
and Tvib.ex. > 2000 K in the disk of HD141569. Note 
that, while the vibrational temperatures of the v=2-l 
transitions are higher than the rotational temperature 
in the UV regime, they are typically much lower than 
the vibrational temperature of even higher-lying levels 
(u=6-5 can be as high as 6000 K). Therefore, while we 
use the v=2-l temperature since it is easy to measure in 
most disks, our results are consistent with previous work 
citing very high temperatures of high-lying vibrational 
levels. 


4.3. An empirical sequence of disk dispersal 

In Figure [6] we illustrate how the temperature-radius 
diagram of rovibrational CO emission resembles key as¬ 
pects of the current theoretical picture of protoplane¬ 
tary disk dispersal. In the phenomenological class of the 
double-component disks, the BC is observed from a hot 
disk region extendi ng to the magnetosp heric accretion ra¬ 


dius at « 0.05 AU (Bouvier et al. 2007), demonstrating 
the presence of abundant gas at the smallest radii dy¬ 
namically permissible for Keplerian orbits. These disks 
are likely to be in a primordial phase, before giant plan¬ 
ets have fully formed (sketch 1). DR Tau is an example 
of a primordial disk as based on its CO emission. 

The next stage in the CO sequence leads to the sup¬ 
pression and eventual disappearence of BC (sketch 2), 
strongly suggesting the opening of partly devoid inner 
gaps or an inner region depleted in gas compared to the 
primordial disks. In fact, it has been observed by moni¬ 
toring CO spectra of a strong episodic accretor, EX Lup, 
that the BC gets significantly weaker relative to NC af¬ 
ter an accretion outburst that depleted th e gas mass in 
the inner disk by an order of magnitude (Banzatti et al. 
2015|). In this stage, double- and single-component disks 


overlap in the diagram. NC is detected from fractions of 
an AU up to « 3AU. At similar radii, and even larger, 
CO is observed in single-component disks, where the BC 
is not detected. TW Hya is an example of this second 
class of disks, which probably has significant overlap with 
the “pre-tr a nsitional” disks as defined from dust emission 


(Espaillat et al. 2007). Indeed, for TW Hya, an inne r gap 


vet et al. 

2002) 

does not c 

lirectl 


Note that a non-detection of the BC 


gas at <0.1 AU. Residual gas and dust in the innermost 
region, which absorbs most of the direct UV radiation 
from the star, may likely explain why the temperature 
profile in some single-component disks continues along 
an IR-pumped power law. Photo-evaporation probably 
plays a role in opening disk gaps, and it is interesting 
to note that UV photo-evaporation is expecte d to take 
over at a few AU (Alexander et al.||2006 2014), close to 
where the transition between the IR and U V regimes is 
found in the CO diagram. Other clearing mechanisms 
may dominate the disk dissipation at smaller disk radii, 
including disk accretion and dust evolution to growin g 
planetesimals and protoplanets (e.g. Gorti et al. 2015). 

The transition to a UV-dominated regime - oF CO rovi¬ 
brational excitati on d efines a third class of disks, as dis¬ 
cussed in Section 32J In order for UV radiation to ex¬ 
cite the gas to high vibrational temperatures beyond a 
few AU, the inner disk region must be largely devoid of 
shielding dust. We identify these disks as having large 
inner gaps, with very little amounts of residual material 
in the inner disk (sketch 3) . This interpretation is c onsis- 
tent with that proposed by | van der Plas et al. (2015), who 
also suggested that UV-pumped CO emissi on may be a 
distinct signature of devoid inner gaps, and Maaskant et 


al. (2013), who imaged large inner dust gaps in two Her 
big disks in the UV branch of our diagram, one of which 
was previously unknown (in HD97048). In total, large 
inner gaps in the dusty disks have been detected in up to 
6 out of 8 disks in the UV branch so far, and for the re¬ 
maining two (HD98922 and HD 190073) the current lack 
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Figure 6. Left: The CO temperature-radius diagram is illustrated in the context of disk dispersal and planet formation as discussed 
in Section^] The data points are identical to those shown in Figure |5l but the stellar masses are indicated by the symbol size. Disks 
with previous detection of gaps in the dust emission are marked with mled symbols. An upper limit on the vibrational temperature for 
rovibrational CO detected in the debris disk aroun d Beta P ic is also incl uded (|Troutman et al. 2011). Right: A sketch illustrating the 
evolution of a protoplanetary disk (based in part on Alexander et al.|2 014), matched to the different regions of the CO temperature-radius 


diagram. 

of resolved imaging may simply be due to their large dis¬ 
tance (> 300 pc). Potentially, our analysis suggests that 
all disks that show a single rovibrational CO component 
are developing or have developed gaps in their inner 10 
AU region, making them prime targets for future imaging 
detections. 

A stellar-mass dependence of Rco may be implied by 
Figure [6j where higher-mass stars tend to have larger 
Rco relative to lower-mass stars. CO gas from < 2 AU 
is rarely detected in disks around stars of > 1.5M©. This 
cannot be explained by stellar magnetospheric truncation 
of the disk, which is expected to happ_en at radii < 0.1 AU 
across the stellar mass spectrum (Bouvier et al. 2007), 
and an inner disk radius set by dust sublimation may 
be marginally consistent only in a few disks, and only 
at < 1 AU given their bolometric luminosities ( Dulle- 


mond & Monnien 2010 van der Plas et al. 2015 k Al¬ 


ternatively, a stellar-mass dependent inner disk lifetime 
would be consistent with the CO temperature-radius dia¬ 
gram, with disks around intermediate-mass stars dispers¬ 
ing faster than those around low- and solar-mass stars 


hs- 

Muzerolle et al.|2010; 

Yasui et al. 

|2014 Ribas et al. 

2015 

). 'That is, we detect i 

ewer broad ( 

IX) components in 


mass-dependence in Figure |H| would be directly linked to 
the timescale for gas dispersal in inner disks, and there¬ 
fore to planet-formation processes and timescales. 

In summary, we suggest that the CO temperature- 
radius diagram shows a sequence of gap opening in pro¬ 
toplanetary disks, where formation of smaller gaps (in¬ 
terpreted as an early transition, in evolutionary terms) 
is marked by the disappearance of the CO broad com¬ 
ponent, while opening of larger gaps (or a more evolved 
phase) is marked by an inversion in the trend of the tem¬ 
perature profile. The UV pumping regime may represent 
the very last step before the primordial gas in the disk is 


dispersed, setting the beginning of the debris disks phase 
(see the Beta Pic point in F igure [6} from rovibrati onal 
CO v=l-0 lines detected by | Troutman et al.||2011|). A 
high-resolution spectrum of CO emission at 4.7 ft m may 
therefore be a powerful discriminator between primor¬ 
dial, pre-transitional and transitional disk stages at radii 
between 0.1-10 AU. Observations of the relative radial 
profiles of CO gas and of dust grains of different sizes 
in inner disks may provide a powerful tool to distinguish 
between competing gap formation mechanisms. 

4.4. Inner disk radii and the distribution of exoplanets 

The lifetime of the planet-forming region of protoplan¬ 
etary disks is linked to the observed orbital distribution 
of exoplanets. In particular, the semi-major axis distri¬ 
bution of giant exoplanets is known to contain structure 
indicating the presence of multiple populations, includ¬ 
ing the “hot Jupiters”, orbi ting very close t o the star at a 
small fraction of an AU (Marc y et al. 2005) , and a d eficit 
of giant planets out to 0.5 AU (Cumming et al. 2008). It 
is generally accepted that hot Jupiters migrated inwards 
from a natal orbit at a few AU, either due t o angular 
momentum exchange with protoplanetary gas (|Lin et al. 
1996 Kley & Nelson 2012) during the lifetime of the disk, 
or due to a later dynam i cal process, such as the Kozai 


1962 Eggleton & Kiseleva-Eggleton 


mechanism (|Kozai 
2001 ) 

To explore this in the context of our observed in¬ 
ner disk radii, we show the distribution of Rin_ com¬ 
pared to the semi-major axes of known exoplanet^Jwith 
Msin(i) > 0.5Mj up in Figure [7] The BC radius distri¬ 
bution matches the orbital distribution of hot Jupiters 
(0.01-0.1 AU), as expected in the Type II migration sce¬ 
nario. That is, we observe that gas is indeed often avail- 


2 From exoplanets.org ( Han et al.|2014| ); 


; June 18 2015. 
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Figure 7. Top: histograms of the inner radii of rovibrational CO 
emission (separated into BC and NC+SC) and of the distribution 
of semi-major axes of known exoplanets with Msin(z) > 0.5 Mj up 
(from the database of exoplanets.org). The exoplanet histogram is 
arbitrarily scaled for comparison to the CO histogram. Bottom: 
same as plot above, but separated into two stellar mass bins. Here 
we show the innermost CO gas in each disk, therefore excluding 
the NC in double-component disks. 
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able at the right disk radii within 1-2 Myr to allow hot 
Jupiters to migrate inwards to the inner gas disk trun¬ 
cation radius. However, hot Jupiters pro bably co mprise 
a small fraction of giant exoplanets (Marcy et al. 2005). 
If the inner disk dissipates rapidly, most planets may 
not have time to migrate, and the orbital distribution 
of “cold Jupiters” will be determined by the inner gas 
radius of transitional disks, as traced by the narrow CO 
component. Indeed, Figure [7] shows that the NC distribu¬ 
tion peaks slightly inward of the bulk of the exo-Jupiter 
distribution (1-10 AU). 

A stellar-mass dependence of the inner disk gas ra¬ 
dius may also be reflected in the exoplanet population. 
It is known that intermediate-mass stars ha ve ver y few 


ho t Jupiters compared to lower-mass stars (Johnson et 


al. 2007 Sato et al. 2008), although a possible overes- 


timation of the intermediat e-mas s stellar masses is cur ¬ 
rently under debate (Lloyd 2011 Johnson et al. 2013). 


Currie (2009) proposed that the paucity of hot Jupiters 
around intermediate-mass stars may be due to a stellar- 
mass dependent disk lifetime, where inner disks around 
intermediate-mass stars dissipate before any giant plan¬ 
ets have time to migrate to the disk corotation radius. 
In the lower panel of Figure [7| we compare the distri¬ 
butions of innermost gas disk radii separated into two 
stellar mass bins at M* = 1.5 M 0 . This distribution 
suggests that the BC component is rarely found in disks 
around intermediate-mass stars, relative to disks around 
low- and solar-mass disks. That is, the paucity of the BC 
within 0.1 AU in disks around intermediate-mass stars 
is reflected by a similar paucity of hot Jupiters in the 


same stellar mass bin. A likely interpretation of our 
CO observations is that the gaseous disk at 0.01-0.5 AU 
around stars with M* > 1.5 M© indeed dissipates faster 
than the planet migration time scale. Larger samples of 
CO rovibrational spectra of disks and exoplanets around 
intermediate-mass stars are needed to confirm this link 
between the natal disk properties and planetary archi¬ 
tectures. 


5. SUMMARY & CONCLUSIONS 

We presented a new analysis of VLT-CRIRES high- 
resolution (R^96,000) 4.7 /im spectra of rovibrational 
CO emission from protoplanet ary disks. We applied an 
empirical decomposition of CO line profiles to the sam¬ 
ple of CO-emitting disks, finding that rovibrational line 
profiles consist of two distinct velocity components, a 
broad (BC) and a narrow (NC) component, primarily in 
disks around solar-mass stars, and of a single CO com¬ 
ponent in transitional and Herbig disks. When splitting 
the observed line profiles into two components, we find 
an empirical relation between the CO line width and the 
disk inclination as expected when the gas is in Keple- 
rian rotation in a disk. This allowed us to estimate the 
inclination for 15 disks in the sample for which it was 
previously unknown. 

Using the profile decomposition to separate blended 
information from different disk radii, we define an em¬ 
pirical diagram of the Keplerian CO emitting radius and 
the vibrational excitation temperature, leading to the 
definition of a universal excitation temperature profile. 
Between 0.05 and 3AU, we find a power law T v ib. e x. oc 
r- o.3±o.i. Beyond 3AU, the vibrational temperature 
turns around and begins to increase with radius. We in¬ 
terpret this as the presence of two different excitation 
regimes, with pumping by infrared photons dominat¬ 
ing in the inner disk, and fluorescence excitation by UV 
photons beyond 3AU. We discuss the CO temperature- 
radius diagram in the context of the theoretical picture 
of inside-out protoplanetary disk evolution, where the 
innermost disk is dissipated first. We find evidence that 
the CO rovibrational line profiles are tracers of the evo¬ 
lutionary stage of inner protoplanetary disks, potentially 
providing a gaseous equivalent to the traditional defini¬ 
tion of classical, pre-transitional, and transitional disks 
as based on the dust continuum emission. Specifically, 
we propose that disks exhibiting only a single component 
in rovibrational CO should be in more advanced stages of 
gap opening and inner disk dissipation than the double¬ 
component disks, where the disappearance of BC marks 
an earlier stage and the UV pumping regime a final stage 
before the debris disks phase. 

Rovibrational CO emission is sensitive to the inner¬ 
most location of abundant gas, including for transitional 
disks, which may have more than one inner region due 
to the presence of gaps. Consequently, we compared the 
observed distribution of inner disk molecular gas to the 
semi-major axis distribution of giant exoplanets. We find 
that the distribution of CO emitting radii for the broad 
component matches that of the semi-major axis distribu¬ 
tion of hot Jupiters. This is consistent with the majority 
of hot Jupiters migrating inward due to Type II migra¬ 
tion, and being halted by the inner truncation of the 
gaseous disk during an early evolutionary phase. The 
radial distribution of the narrow CO component peaks 
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at larger radii, close to the “period deficit” between the 
hot Jupiter population and the bulk of the giant exo¬ 
planet distribution. Finally, we find that disks around 
intermediate-mass stars rarely have gas in the hot Jupiter 
region. This matches an observed deficit of hot Jupiters 
around intermediate-mass stars, lending support to the 
idea that their inner disks may evolve too fast for giant 
planets to migrate inwards to the disk corotation radius. 

This analysis ultimately supports the potential of fu¬ 
ture spatially-resolved imaging of rovibrational CO to 
move significant steps forward in our understanding 
of disk evolution and planet formation in protoplane¬ 
tary disks. Herbig disks can already be resolved with 
adaptive-optics assisted observations, for instance with 
CRIRES, but disks around solar-mass stars require ei¬ 
ther spectro-astrometric techniques, infrared interfero¬ 


metric imaging (e.g., VLT-MATISSE), or full aperture 
spectro-imaging at 4.7 fi m with the future generation of 
extremely large telescopes (e.g., E-ELT-METIS). 
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tary Geology and Geophysics Program under grant NAG 
5-10201. This work is based on observations made with 
ESO telescopes at the Paranal Observatory under pro¬ 
grams 179.C-0151 and 093.C-0432. 


APPENDIX 


CO VIBRATIONAL RATIOS 

The vibrational ratio v2/vl is the flux ratio F v= 2 -i/F v= i_q between the v = 2 — 1 lines and the v = 1 — 0 lines 
with same rotational quantum numbers in the two vibrational levels. In this work, v2/vl values are measured from 
stacked lines built by using multiple rotational lines and thus provide an average estim ate of the v ibr atio nal ratio 
between v = 2 and v = 1. Previously, vibrational ratios of CO emission were measured byjBrown et al.| (2013) for 25 
disks that are included in our sample, using the same CRIRES spectra we use here. They published the median of 
five values measured in each disk, as obtained from ratios of the co rresp onding rotational lines in the two vibrational 
levels. Figure |8| (panel a) compares the values from Brown et al.f (|2013|) to those measured in this work. The values 
measured in the single-component disks agree to a 1:1 relation, showing that our estimates are consistent. However, 
the vibrational ratios in previous work have been underestimated in the double-component disks, because the different 
contributions from BC and NC in the v=l-0 lines were not separated. 

The vibrational ratio v2/vl depends on the temperature and column density of the gas, as well as on the excitation 
process. In this work we use a model of a slab of gas in thermal equilibrium that accounts for the line op acity, where 
the level populations are determined by the Boltzmann distribution (the model is explicitly described in Banzatti et 
ah| 2012). The mo lecular p arameters for CO rovibrational lines are taken from the 2012 HITRAN and the H1TE1V1P 
databases (Rothman et al. 2010,2013). We produced a grid of models to convert the measured v2/vl into an estimate of 
the vibrational excitation temperature T v ib. ex . of the gas. We used three values for the CO column density N co t o span 
the range of values proposed by previous studies of NIR CO emission from protoplanetary disks (10 17 -10 19 cm 


Blake 


& Boogert|2004| |Herczeg et al.|2011| |Salyk et al.|2011l|Banzatti et al.|2015| |van der Plas et al.|2015|) . Similarly, tor the 
temperature we use eleven values between 300 and 5000 K. Lor each set of 'iVib.ex. and Aqo a CO model spectrum is 
produced, assuming a line FWHM of ~ 75 kms -1 (the median value for the BC) and using a pixel sampling equivalent 
to the CRIRES spectra. We then re-sample the model spectrum using a typical value for the rms noise as measured 
in the observed spectra, and then measure vibrational ratios following three methods: 1) using directly the individual 
line fluxes before convolution with the observed FWHM (the “true” value), 2) taking the median of seven ratios of the 
line fluxes measured in the convolved spectrum (equivalent to the method used by Brown et al. 2013), and 3) using the 
method adopted in this work. Figure M (panel b) shows one model as an example, for illustration of the outcomes of 
this procedure. It is important to highlight that the “true” ratio is not a single value. In the example, the vibrational 
ratio varies between 0.2 and 0.6 depending on the line used to measure it. Moreover, vibrational ratios measured from 
the individual line fluxes as observed in the convolved spectrum tend to underestimate the true value. This is due 
to blending of the v=l-0 lines with higher vibrational lines, which contribute increasingly to the line flux as T v ib. e x. 5 
Nco, and the line FWHM increase. The averaged ratios estimated in this work in general agree well with the range 
of true values for each model (panel c in Figure [8]). This is true in particula r in th e case of a CO column density of 
10 18 cm -2 , which is a good approximation for most disks in our sample (e.g. S alyk et al.| 2011). For higher values of 
Nco = 10 19 cm -2 , the vibrational ratios measured here may overestimate the true values by up to 50%. 

We convert the vibrational ratios measured in this work into temperatures using the curves shown in panel c of 
Figure [8j In this regard, it is not important that our measured v2/vl may overestimate the true value, as long as 
we are self-consistent in the conversion, i.e. as long as we use the same method to measure v2/vl in the data and 
in the model. In Tables [ 2 ] and [ 3 ] we put the temperature values estimated by using a CO column of 10 18 cm -2 . In 
the case of AS 209, HH 100, and IM Lup, the measured vibrational ratios are too high to be modeled with a CO 
column of 10 18 cm -2 , which never produces values higher than ^0.8. Therefore, for these disks we use a CO column of 


10 19 cm -2 , which can reproduce v2/vl as high as 1.2 and would be appropriate if the emission is more optically thick 
(our analysis indeed finds that the BC in AS 209 and IM Lup is more optically thick than in other double-component 
disks, as the BC is detected also in 13 CO, see Figure [ 2 ]). For each one of the three values adopted for the CO column 
density, we obtain a different empirical relation between the temperature and the radius of the CO emission. This is 
shown in panel d of Figure [8j where we plot the linear fits obtained as in Figure [4j The index q of the temperature 
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Vibrational ratios T = 1000 K , logN = 18 cm 2 






Figure 8. In order from top left to bottom right. Panel a: comparison of vibrational ratios measured in this work to those measured in 
|Brown et" al. ( |2013| ). The 1:1 relation is marked with a solid line. Panel b : comparison of measured values for v2/vl in a model with T = 
1000 K and N = 10 18 cm -2 , showing the range of P-b ranch rovibrational levels covered by the CRIRES spectra; the true value is shown 
with a solid black curve, while the method adopted in|Brown et al.|(|2013|) and in this work are shown in grey-dashed and red-solid lines 
respectively. Panel c : vibrational ratio as a function of temperature and column density over the grid of models explored; the range of true 
values for each model is shown as a vertical thick solid line. Panel d: radial temperature profiles obtained after conversion of the measured 
v2/vl values into temperature values as shown in panel c for three values of the CO column density. An upper limit to the index q is 
overplotted in dashed-dotted line, consistent with the data in the case of a gradient of CO column densities increasing with disk radius. 


profile T v ib.ex. oc R q varies between 0.20 (for Nqo = 10 19 cm 2 ) to 0.33 (for Nqo = 10 17 cm 2 ). If we allow a radial 
gradient in TVco such to have higher values of ~ 10 18 — 10 19 cm -2 as representa tive of the m o re opt icallv-thick NC, 
and lower values of ~ 10 17 — 10 18 cm -2 as representative of the BC (as found by Herczeg et al.||2011), the index may 
become as steep as rsj 0.5. 


CO DECOMPOSITION OF THE ENTIRE DISK SAMPLE 

In Figure [9] we show the outcome of the CO decomposition procedure as applied to all disk spectra in our sample. 


BAYESIAN METHOD FOR LINEAR REGRESSIONS 


We i mplement the Bayes ian method by Kelly (2007| in the linear fits to the empirical relations presented in Sections 


3.3 and 3.5 Kelly (2007) provided a Bayesian method for linear regression that accounts for errors on both the 


dependent and independent variables as well as an intrinsic scatter due to physical properties that are not explicitly 
included. A weighted sum of Gaussian distributions is used to obtain an accurate approximation of the true probability 
distribution of the independent variable. To sample the posterior distribution we use the Markov chain Monte Carlo 
method from the Metropolis-Hastings algorithm, which is more efficient than the Gibbs sampler in case of small s ample 
sizes. A thorough discussion of the method and of its comparison to other statistical methods is provided in |Kelly| 
([2007). The method has been made publicly available through the IDL routine linmix_err .pro, that we use in our 
analysis. Figure [lO] shows the posterior distributions we obtained for the linear regression parameters presented in 
Section [3j 
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Figure 9. CO decomposition plots of the entire disk sample. The color-coding follows what used in Figure [ 2 ] 12 CO v =1—0 lines are 
shown in black, 12 CO v = 2 — 1 lines in orange; NC is filled with cyan horizontal lines, BC with red diagonal lines, SC with blue horizontal 
lines. In the double-component disks (first three columns from the left), 12 CO v = 2 — 1 lines are scaled up to the 12 CO v = 1 — 0 lines to 
illustrate the two CO components. In the three disks where NC is detected also in the 12 CO v m 2 — 1 - DRTau, AS205 N, and SCrA N - 
Gaussian fits to BC and NC are shown on top of the line. Single-component disks are shown in the last two columns to the right. 
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Figure 10. Posterior distributions for linear regression parameters obtained with the method by [Kelly (2007| for the empirical relations 
presented in Section [3] The first two plots from the top are for the empirical relations between the BC and NC FWHM and the disk 
inclination (Figure [ 3 ]), while the last two plots are for the relations between vibrational ratios, or temperature, and disk radius (Figure [4|. 
Vertical lines mark The median value of each distribution, which is reported inside each plot together with a robust estimate of the posterior 
standard deviation. 
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